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Abstract 

The effect of a sinusoidal modulation of the interaction strength on a Fermion pair condensate 
is analytically studied. The system is described by a generalization of the coupled fermion-boson 
model that incorporates a time-dependent intermode coupling induced via a magnetic Feshbach 
resonance. Nontrivial effects are shown to emerge depending on the relative magnitude of the 
modulation period and the relaxation time of the condensate. Specifically, a nonadiabatic modu- 
lation drives the system out of thermal equilibrium: the external field induces a variation of the 
quasiparticle energies, and, in turn, a disequilibrium of the associated populations. The subsequent 
relaxation process is studied and an analytical description of the gap dynamics is obtained. Re- 
cent experimental findings are explained: the delay observed in the response to the applied field is 
understood as a temperature effect linked to the condensate relaxation time. 
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The study of ultracold atomic gases has led to a remarkable series of experimental real- 
izations of fundamental effects Essential to many of these achievements has been the 
control of the interaction strength via a Feshbach resonance (FR), which has allowed the 
emergence of these systems as a practical testing ground for quantum-statistical and many- 
body physics. Specially relevant has been the realization, with a two-component Fermi gas 
of atoms, of the crossover from a molecular Bose-Einstein condensate (BEC) to a Bardeen- 
Cooper-Schrieffer (BCS) superfluid of loosely bound atom pairs 

By applying a magnetic FR, the interaction is led to change from repulsive in the BEC phase 
to attractive in the BCS side. Whereas two-body physics supports a bound molecular state 
in the BEC side, the formation of pairs in the BCS regime occurs only due to many-body 
effects. Significant advances have been made in the understanding of different aspects of 
;his transition. In particular, the role of thermal fluctuations has been extensively analyzed 



12l |. Despite those advances, further work on the characterization of nonequilibrium aspects 



of the crossover is required. Indeed, the extension of the experiments on scattering-length 
variations to unexplored time-dependent regimes and the setup of an expanded theoretical 
framework where the emergent effects can be understood are current challenges. Here, we 
aim at explaining recent experimental results on the nontrivial dynamics resulting from a 
sinusoidal modulation of the interaction strength [l^. This will require a precise character- 
ization of the different time scales, in particular, of the relaxation time of the condensate. 
The interest of the study is not restricted to the field of ultracold atomic gases. In fact, 
the expansion of the BCS model proposed in our approach can have broad applicability: 
the effect of a disequilibrium population on the gap dynamics, which is the central issue in 
the analysis, is relevant to topics ranging from nonequilibrium superconductivity (l^ . [isl to 



16|. 



quench dynamics in superfluid ^He 

As starting point we take the experiments reported in Ref. [13]. In them, the interaction 
strength of a gas of ultracold ^Li atoms in the BCS regime was sinusoidally modulated 
through a magnetic FR. (The (broad) FR at 834 G between the two lowest hyperflne states 
was used.) The system response, which consisted in a damped oscillation of the condensate 
fraction with the modulation frequency, was found to be delayed with respect to the applied 
field. The delay showed no appreciable changes at different cycles of the external field; 
moreover, it presented the same scale for widely different frequencies. Additionally, the 
damping time was observed to be much longer than the driving period. In a preliminary 
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analysis, the deferred response was conjectured to be rooted in the finite relaxation time of 
the condensate; furthermore, the decay of the oscillation amplitude was linked to heating 
resulting from the nonadiabaticity (on the gap time scale) of the process. (In the present 
paper, "adiabatic" will also be applied to a modulation much slower than the condensate 
relaxation.) To evaluate those conjectures, the measured characteristic times were compared 
with related theoretical predictions. However, that analysis was not conclusive about the 
origin of the recorded behavior because of the limitations of the available models, specially, 
of the lack of an appropriate description of finite-temperature effects. Our objective is 
to provide a theoretical framework where the observed features can be understood, and, 
in particular, the previous conjectures can be assessed. To this end, we concentrate on 
conditions which allow an analytical description of the dynamics, and, consequently, a clear 
identification of the dominant mechanisms. 

We consider a gas of ultracold Fermi atoms with two hyperfine states coupled to a molec- 
ular two-particle state via a magnetic FR. Our methodology to deal with magnetic-field 
modulations combines three main elements. First, the standard framework, namely, the cou- 
pled fermion-boson model 17, Q, 1^, [2^, is expanded by incorporating a time-dependent 
intermode coupling. Second, a perturbative scheme, valid for a sufficiently small modu- 
lation amplitude, is introduced in the Hartree-Fock-Bogoliubov (HFB) description. In this 
approach, finite-temperature effects are tackled and the out-of-equilibrium situation induced 
by the modulation is characterized. Finally, a method for describing the evolution of the 
quasiparticle populations, and, subsequently, the gap dynamics is presented. Accordingly, 
we start from the unmodulated system: the grand-canonical Hamiltonian reads 



H - fiN = gkQk.gQk.g + Vint ati_^^^ati_,^^^aa_k',iaf +k',T + 

k,CT q,k,k' 

J2 (^q + ^^0) &i&q + gYl (^q«H+k,T«a-k,i + h-c.) (1) 

q q,k 

where /i is the chemical potential, is the total number of bare Fermi atoms, aj^^ (ctk.cr) 
denotes a fermionic creation (annihilation) operator of an atom with momentum k and spin 
cr, (cr G {Tj i})? and 6]^ (6q) is a bosonic operator that creates (destroys) a molecule with 
momentum q. The populations corresponding to the two hyperfine states are assumed to 
be equal. The free dispersion relations for fermions and bosons are = /?.^k^/2m — 
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and = h^q^ /Am — 2/i, respectively. The binary attractive interaction potential between 
fermions is characterized by ^nt(<0). Additionally, g represents the FR coupling between 
the closed and the open channel states, being the detuning of the boson resonance state 
from the collision continuum. 

Initially, the system is at equilibrium at a finite temperature T. In that situation, a 
sinusoidal modulation of the detuning from the FR is applied. Correspondingly, is replaced 
by v{t) = z/q + Asmujpt. It is assumed that Vint^ which characterizes the attractive pairing 
interaction resulting from nonresonant processes, is not affected by the applied detuning of 
the FR. Through the unitary transformation U{t) = e* ^'^^'^ , the Hamiltonian is 

transformed into H' = WHU — ihWU; consequently, Eq. (1), [with z/q replaced by z^(t)], is 
rewritten as 



H' — fiN = e^al^^^ai,^^ + Vint ^ ctq+k,T'^a_k,i^2~'^'>i'^2+i'''f + 
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k,(T q,k,k' 



q V q,k 



Our procedure to analyze the dynamics starts, like the standard HFB approach 
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with the introduction of three mean fields: n = X^k (|'^k cr'^k.a^ for the spin density, 
A = \Vint\ Xlk ('^-k,i'^k,T) for the pairing field, and 0^ = (^q=o) for the boson field. (We 
take q = as we focus on the condensed molecular field.) Next, a perturbative scheme 
is set up as follows. From Eq. (2), it is apparent that the effect of the magnetic modu- 
lation can be understood as a time variation of the intermode coupling strength: we can 
work with the effective strength gefj = ge^p cosujpt _ ^ ^ ^^^^^^ where Sg{t) stands for the 
modulation-induced increment. (Note that the time dependence of geff prevents the effective 
one-channel reduction applicable, for a broad FR, to the undriven coupled fermion-boson 
model.) Furthermore, from the expansion e* = ^'^_^i^JiiA/ujp)e''''^p^, and taking 
into account the properties of the Bessel functions, it follows that, for A/up <^ 1, we can 
make the approximation 6g ~ i2gJi{A/ujp) cosupt, the magnitude of 6g being much smaller 
than that of g. (Higher-order terms will be discussed later on.) In turn, the previously 
defined mean fields can be expressed as n = uq + 6n, A = Aq + SA, and (pm = 4>m,o + S4>m, 
where no, Aq, and 0^,0 are the respective values in the absence of the magnetic variation, 
and, 6n, 5A, and 6(j)m are the corresponding modulation-induced increments. To first order. 



the complete Hamiltonian can be split as H' — fiN ~ Ho + Hper. The zero-order term, which 
describes the unmodulated system, is given by Hq = a ^^tjc cr'^k.cr — Xlkl^octk |Ct^k j,+h.c.) 
and corresponds to an effective BCS model with mode energy V^. = — + Vintno and gap 
Ao = Ao — g4>m,o- The first-order correction reads Hper = —6A{t) J^k '^k t'^Lk j, + h-c, where 
6A(t) = 6A—g6(l)m—Sg(t>mfl is the increment of the generalized order parameter A = A— gf^m- 
(We have neglected the variation in the atomic density, i.e., we have taken 6n ^ 0, which 
is justified for a broad FR 0|. The generalization required to deal with a narrow FR is 
straightforward.) To complete the mean-field description, we must add the equation for the 
evolution of the boson mode, namely, ifi^^^^^^^^^ = {1^0- 2^) (0„,o + + ^ ( Aq + 5 A) , 
which, through the application of our perturbative approach and taking into account that 
the equilibrium molecular field is given by 0m,o = |v- tf(tu-i^o) ' converted into 



= (^0 - 2/i)50„ + -^5A + Ao. (3) 

"J I yint I I yint \ 

Here, the presence of the driving term p^I^Aq points to the oscillation of 50m with frequency 
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Hq is standardly diagonalized through a Bogoliubov transformation (BT) characterized 
by the fermionic operators Ck,t = cos6'fcak,T ~ sin^^^j^^ and clj^^ = sin^fcCk,! + cos6'fca^,^|, 
where Ok is defined by tan(26'fc) = Aq /Vk Q, With this BT, Hq is cast into = 
Xlk -^fc.o(Ck |Ck,| + cj^|Ck,|) + constant. Importantly, cj^ ^ (ck.t) corresponds to the creation 
(annihilation) operator of a quasi-particle excitation with momentum k and spin | from 
the BCS state, which acts as an effective vacuum state. The associated excitation energies 
are E^fi = \Jv^ + Ag. The excitation gap Aq combines the mean field Aq, obtained from 
the BCS equation Aq = X]k(2/fe " 1) sin26'fe, with the equilibrium molecular field 0^,0, 
which results from the equation for the boson mode in the absence of driving. As thermal 
equilibrium is assumed for the system before the application of the magnetic modulation, 
the populations of the quasiparticle states {fk} are initially given by the Fermi distribution 
function fll = 1/(1 + e^^^«/'=s^). 

Through the previously defined BT, the perturbation Hamiltonian is converted into 



Hper = ^SEk{t) (^4^^Ck,t + Ck,i^k,i) + k{t)cl^^cl^^^ + h.C.j , (4) 
k 
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where SEk{t) = i(5A(t) sm26'fc + c.c, and Gk{t) = -(5A(t) cos^ 6*^. + 6A*{t) sin^ Ok. From 
the form of Hper, a first picture of the dynamical implications of the field modulation can 
be drawn. The (time-dependent) diagonal terms lead to a time variation of the quasi- 
particle energies, which become Ek{t) = Ekfi + 6Ek(t). The non-diagonal terms represent 
modulation-induced interactions between the vacuum state and a doubly-excited state. Im- 
portantly, these coupling terms, which oscillate with the external frequency oOp, are relevant 
only when they can induce an effective resonance between the BCS state and the two- 
excitation configuration, i.e., only when ujp > 2Ao, {h = 1.) Here, in order to isolate the 
mechanisms responsible for the delayed response, we concentrate first on the regime defined 
by Up < 2Ao. In this frequency range, the interaction terms can be discarded, and the per- 
turbation Hamiltonian can be approximated as -f/per = X^k '^-^^('^k t'^k,? + '^kj'^k.i)- Hence, 
the quasiparticle states of the unmodulated system still provide a diagonal representation 
of the driven Hamiltonian. The appearance of heating effects outside this regime will be 
discussed later on. 

Now we turn to incorporate finite-temperature effects in the above framework. The 
modulation drives the system out of equilibrium as the initial populations, i.e., the thermal 
values associated with the unmodulated energies, do not correspond to the Fermi distribution 
/fc^(^) ~ V(l +6^'=*^*-'/'^^^) for the actual (time- varying) energies. The description of the 
effect of this quasi-particle disequilibrium on the gap dynamics requires a selfconsistent 
approach since the energies and the gap are interdependent. Indeed, as shown by the 
expression Ek(t) = E^"^ + {^-;^^5K{t) + c.c. j , the quasi-particle energies are affected by the 
gap evolution an d by the molecular-field variation; in turn, the {Ek(t)} enter the general 
gap equation 2^, 



A(t) =\^Y. [2^W - 1] sin2^fc, (5) 

k 

via the (changing) associated populations {fkif)}- We will see that it is precisely the evolu- 
tion of the populations, more specifically, their relaxation towards equilibrium, that gives the 
keys to understanding the experimental results. An important aspect of this problem can be 
understood by now: a finite relaxation time r/ of the {fk} is necessary for the appearance of 
the gap delay. In fact, for a sudden relaxation, the populations follow adiabatically (on the 
relaxation time scale) the equilibrium values {fl'^if)} associated with the time-dependent 
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energies. The evolution corresponds then to a sequence of equilibrium states where time 
enters as a parameter, the associated gap dynamics being "trivial": no delay between the 
gap evolution and the external field emerges. Therefore, to reproduce the delayed response, 
we must go beyond that adiabatic regime. Accordingly, we present a self-contained deriva- 
tion of the dynamics with no constraints on time scales. The evolution of the populations 
is assumed to be governed by the equation [l^ 



where 1/r/ represents the effective thermalization rate. The relaxation mechanism can be 
conjectured to be rooted in collisions between excited particles. Here, we do not go into 
details of the dependence of Tf on the system characteristics; instead, as we focus on general 
aspects of the role of the population thermalization in the condensate relaxation, we consider 
a generic Tf. Here, a comment on the stability of the temperature is in order. One must 
take into account that the standard trapping conditions allow assuming the stability of 
the temperature for the considered small variations of the scattering length. In this sense, 
we recall that a grand-canonical description, which is routinely applied in this context as 
it can incorporate the possible exchange of particles between the condensate fraction and 
the thermal cloud, implies that a fixed temperature can be reasonably assumed. Also, 
it is worth stressing that a parallel treatment of the fiuctuations of the condensate field 
would be necessary to formally complete our description. However, it is shown that, in 
our perturbative regime, the non-condensate fraction has a second-order effect on the gap 
dynamics. (For a systematic treatment of different aspects of the role of fluctuations, see 



Eq. (6) is exactly solved to give fk{t) = fl\t) - j[^e-'^*-*'^l^f^{t')dt'. Next, this 
expression for the populations is introduced into Eq. (5) to give the following integral- 
differential equation for the order parameter 



dfk 
dt 



(6) 



Ref. [12].) 



Ao + 5A(t) =^ J] 2 




dt 



■eq 



{t')dt' 



) 



1 sin20fc. 



(7) 



Here, one must take into account that f^'^{t) contains 6A(t) and 50m, and contains SA 
and S(j)m- Hence, we have obtained a description of the gap evolution, albeit in implicit form. 
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Given the complexity of this picture, the problem of identifying the origin of the deferred 
response, reduced at this point of the study to that of uncovering the connection between 
the delay time and r/, is still nontrivial. However, it simplifies considerably in the following 
regime, where an explicit characterization of the gap evolution is feasible. Specifically, for 
Ek ^ A <^ T ^ Tc, {kB = 1), where Tc is the temperature for the BCS transition, we 
can make the approximations f^'^it) ~ /^g + j^5Ek{t) and ~ Then, with the 

expression for the unperturbed gap Aq = ^^^^ ^j^. [2/^ ^ — l] sin 29^ and the approximation 
^ ~ -^^f - = sin24((5A - g54)^), Eq. (7) is cast into 

5A = ^sin2 2^fc {6A-g6(Pm)- [ e'^'-'">/^f {6A - g6<Pm)dt' . (8) 

Now, following a standard procedure, this integral-differential equation is converted into the 
differential equation 

—TT = 9Hm, (9) 

dt Tf Tf 

where x = ^'^^W' lo^ ^^^^ 26kk'^dk encapsulates the overall effect of the quasi-particle states 
on the gap response. {K is the upper limit of the momentum summation required by the 
standard renormalization procedure.) In the considered regime, x — ^^^^jP- <^ 1 [22||, where 
hp is the Fermi wave number and a is the background scattering length. Eq. (9) along with 
Eq. (3) for the molecular field constitute a closed set of equations for the system evolution. 
From them, it is apparent that 6A is determined by the combination of effective driving, 
coming from the term ^gdcf)^, and damping with rate As the driving is continuously 
taking the system out of equilibrium, the relaxation mechanism is permanently activated. 
The combined effect of both mechanisms can be expected to produce a nondirect following 
to the external field. Approximate analytical solutions confirm these predictions: the gap 
evolution is given by 

A(t) 



A 







1 -C [e~*/^« sin v? + sm{upt - ^)] , (10) 



where C = 2xtS—! — ^lMZ^£L=(^> g) determines the amplitude of the induced oscillations, 
(p = aictan^LjpTf) is a phase shift with respect to the applied magnetic field, and tr = 
appears as the condensate relaxation time. (The meaning of r^j becomes evident in a 
simplified scenario: for the system with no external driving, a sudden perturbation of the 
gap is shown to relax to equilibrium with characteristic time r/j.) In the considered regime, 
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namely, near the critical temperature and for a perturbative gap variation, it is found that 
Tjic:^ Tf. The correspondence of these results with the experimental findings is summarized 
in the following points. 

(i) The system response contains a transitory decay with characteristic time tr-, and, 
as observed in the experiments, a secular oscillatory behavior with ujp. The amplitude, 
which combines in a nontrivial way parameters of the external field and characteristics of 
the unperturbed system, refiects the complex character of the driving mechanism. The 
following to the external field is not instantaneous: there is a delay time associated with the 
phase shift ip and given hy td = ^ = tr[1 + O {{upTrY)]. Hence, as conjectured in Ref. 
[13], td approximately corresponds to the condensate relaxation time. The delay presents 
no changes at different cycles of the external field. Furthermore, the small magnitude of the 
correction O {{ujpTjiY) for the conditions of the experiments explains the detected invariance 
of the delay scale with the modulation frequency. As refiected by the minus sign before C, 
there is an extra phase shift vr between the gap oscillation and the driving field. This 
corresponds exactly to the results presented in Figs. 2 and 4 in Ref. [13], where an inversion 
of the magnetic-field axis was introduced to facilitate the observation of the delay. 

(ii) The connection between the different time scales is uncovered. Since tr gives the 
time for the condensate to reach the thermal equilibrium, it is directly related to the ther- 
malization rate of the populations. Notice that tr can significantly differ from outside 
the considered regime with x ^ 1. This can be understood taking into account the intri- 
cate interdependence of the gap and the populations, which implies, in general, a complex 
nonlinear contribution of the populations to the gap relaxation |14|. 

(iii) The mechanism responsible for the deferred response is rooted in the finite reaction 
time of the gas to a variation in the quasiparticle energies. [Note that the adiabatic limit 
corresponding to a sudden relaxation, i.e., to Tf 0, is consistently recovered in Eq. (10).] 
Moreover, as conjectured in a preliminary analysis, the observed delay is a temperature 
effect. At zero temperature, there is no initial population of the excited states; furthermore, 
as the modulation in the assumed regime does not induce a transfer from the fundamental 
state, the excited states are never populated. Hence, at T = 0, there is no population 
relaxation, and, consequently, no delay in the gap evolution [23|. 



(iv) It is worth discussing the effects that can be expected outside the considered regime 
of nonexciting frequencies {up < 2Aq) and perturbative amplitudes. First, for ujp > 2Ao, the 
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magnetic field can induce an effective resonance between tfie vacuum state and a doubly- 
excited state. As a consequence, the interaction between those states, represented by the 
non-diagonal terms in Eq. (4), becomes important, and significant heating can result. (See 
Refs. [24] and [25] for related work.) Second, as the amplitude is increased, the contribution 
of the terms of order higher than one in the expansion of the exponential in Eq. (2) grows. 
Given that the frequency of each term is a multiple of cUp, a resonance between the ground 
state and the excitations can eventually be reached for increasing order, which, again, can 
lead to an irreversible loss of population. 

At this point some aspects of our approach must be recalled. Importantly, a perturbative 
regime has been considered: the system, which is initially in the BCS side, is assumed to be 
inside that regime during the whole modulation process. Hence, the unitary limit of large 
scattering length is never reached and neither is attained the EEC side. Our approach, 
directly set up from the fundamental theory, has some similarities with former studies in 
superconductivity where the effect of a population disequilibrium on the gap dynamics was 
tackled by introducing an operative changing temperature in the static Ginzburg-Landau 
(GL) equation Q, [l^. The effective time-dependent GL equation thus obtained was shown 
to satisfactorily explain the relaxation process. That description, as ours, is basically built 
from the incorporation of formal solutions for the evolving populations into the gap equation. 
In the uniform case and in the perturbative regime for the gap variation, the analogy with 
our self-contained approach is complete. Incidentally, we stress that the success of our 
uniform description in reproducing the experimental results, which, in fact, were obtained 
for a harmonic confinement, reflects the robustness of the identified physical mechanisms 
against spatial non-uniformities. The generalization through a local-density approximation 
is straightforward. 

In summary, we have presented an analytical explanation for the delayed response of a 
Fermi condensate to a modulation of the interaction strength. Although a more quantitative 
comparison with the experiments of Ref. [13] requires additional information on the temper- 
ature and the amplitude of the applied field, the study uncovers fundamental aspects of the 
out-of-equilibrium dynamics of the condensate, in particular, the nontrivial role of the relax- 
ation in a time regime of experimental and theoretical interest. Our analysis can have direct 
practical implications. One of the main motivations for the experiments was the validation 
of the detection schemes based on a projection of the Fermion pair condensate into a molec- 
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ular condensate. Those schemes are applicable only if the response time of the condensate 
to variations in the interaction strength is much larger than the sweep time. The relevance 
of the nonzero-temperature character of the delay to the projection techniques is clear: an 
out-of-equilibrium situation and the subsequent relaxation process of the condensate can be 
induced not only by changing the temperature but also by manipulating the system with 
external fields. Our picture provides the theoretical basis for the design of methods for 
measuring the different time scales, and, consequently, for defining the appropriate ranges 
for the projection. Furthermore, the identification of the dominant mechanisms opens the 
way to develop variations of the basic arrangement as probing tools for different aspects of 
the dynamics. Given the generality of the applied model, the applicability of the study in 
parallel contexts can be expected. 
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